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Closed expressions are derived for resonant multidimensional X-ray spectroscopy using the quasi- 
particle nonlinear exciton representation of optical response. This formalism is applied to predict 
coherent four wave mixing signals which probe single and two core-hole states. Nonlinear X-ray 
signals are compactly expressed in terms of one- and two- particle Green's functions which can be 
obtained from the solution of Hedin-like equations at the GW level. 



I. INTRODUCTION 

With the advent of femtosecond to attosecond X-ray 
sources [H, 0, H, 0,[1,Q, time resolved resonant core level 
X-ray spectroscopy [7] has become a promising tool for 
studying electron and nuclear motions in real space and 
real time. Resonant techniques provide detailed infor- 
mation beyond the charge density derived from time- 
resolved diffraction. Nonlinear X-ray spectroscopy can 
monitor the dynamics such as making and breaking of 
chemical bonds at the molecular level. 

In a recent work Q| we derived closed expressions 
for attosecond coherent stimulated x-ray Raman signals 
(CXRS). It was shown that CXRS can be expressed in 
terms of one-particle Green's functions which depend 
parametrically on the core-hole. In this paper we ex- 
tend this approach to study nonlinear x-ray spectra by 
combining the many-body Green's function formulation 
with the exciton representation of valence electronic ex- 
citations. The exciton representation has been par- 
ticularly useful for computing the o ptic al valence ex- 
citations of molecular aggregates [9l lid Hlj . Hedin's 
Green's function approach [12j | provides an exact many- 
body formulation which avoids the computation of many- 
body wavefunctions. In most applications, Green's func- 
tion theory requires a self-consistent computation which 
can become quite tedious and demanding as the sys- 
tem size grows. Various approximations are therefore 
invoked in practical applications. The GW approxima- 
tion [l2|, [H[ constitutes a good compromise of accuracy 
and cost, and various variants have often been used in 
many applications of the Green's function approach. By 
combining the exciton modelling of nonlinear response 
with many-body formulation, we can reduce the com- 
plexity of the Green's function approach. We derive 
a Bethe-Salpeter (BS) type equation for two electron- 
hole pairs. In the standard BS equation [l]|, the scatter- 
ing matrix which represents the interaction between an 
electron and a hole, is expressed as a derivative of the 
one-particle self-energy with respect to the one-particle 
Green's function which is then calculated using the GW 
approximation [14l . Il5| . A further simplification is in- 
troduced by taking the screened-coulomb interaction to 
be independent of the Green's function Despite the 



lack of a rigorous justification, this approach works well 
for a broad range of systems including isolated local- 
ized atoms and molecules [l5j], bulk semiconductors and 
insulators [H, [l7l|, conjugated polymers [l8j and infinite 
periodic crystals [3, [l5| . It has been difficult to quan- 
tify the influence of the variation of screened potential 
with respect to Green's function on the optical spectra 
because it involves an expensive self-consistent calcula- 
tion of the derivative of the screened potential. Here we 
employ the GW approximation to derive an algebraic 
expression for the electron-hole scattering matrix which 
only depends on the one-particle Green's function. The 
variation of screened-coulomb potential is included ap- 
proximately. This yields a simplified practical expression 
for the scattering-matrix which requires comparable ef- 
forts to computing the one-particle Green's function at 
the GW level. We extend this result and combine it with 
the nonlinear exciton equation (NEE) approach [l9|,|20(to 
compute nonlinear X-ray signals. 

The paper is organized as follows. In the next section, 
we introduce the one- and two-particle Green's functions 
parametrized in terms of core-holes. By defining a pro- 
jection operator in the one-excitation space, we derive a 
closed expression for the nonlinear signal in the direction 
km = ki+k2— k3 in terms of parametrized Green's func- 
tions. In Sec. piip. we present a self-consistent scheme 
for computing these Green's functions. We conclude in 

Sec. eg. 

II. GREEN'S FUNCTION EXPRESSION FOR 
THE NONLINEAR RESPONSE 

We divide the electronic system of the molecule into 
a core and the valence parts. The system is described 
by the Mahan-Nozieres-De Dominies (MND) deep core 
Hamiltonian [U [H, [1| 

H = H + H c + Hi (1) 
where Hq is the free electron part, 

H = ^2 + e « £ n 5 «- ( 2 ) 

i n 

The indices i,j,k,l denote the valence orbitals whereas 
m, n are the core orbitals. c\ (a) are Fock space Fermi 
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creation (annihilation) operators. The coulomb interac- 
tion is 

iji'f 
ranm'n' 

— WimnjclcmC^Cj . (3) 

imnj 

The three terms represent respectively valence- valence, 
core-core and valence-core interactions. The dipole in- 
teraction with the X-ray pulse, 

Hi = - J2 ( E+{t)^ m c\c m + E- (t^clci ) (4) 

im 

where the field is treated as a scalar for simplicity. We 
assume that the core-holes created by exciting core- 
electrons to the valence orbitals have an infinite mass. 
Their motion can thus be ignored and they enter the 
calculation as fixed parameters. Only valence electron 



dynamics needs to be considered. This is a reasonable 
approximation for X-ray spectroscopy where a hole is 
created in one of the core orbitals which is tightly bound 
to the atom nucleus compared to the valance electrons. 
Core mig ration considerably slows down for deeper core 
levels |26j. 

Using the exciton quasiparticle representation, the co- 
herent nonlinear response functions can be expressed in 
terms of the exciton scattering-matrix which represents 
the interaction between two excitons created b y in terac- 
tions with X-ray pulses at different sites (Til. l3l|. For 
clarity we focus on the double quantum coherence sig- 
nal generated in the direction km — ki + k2 — k3 [1JJ . 
Other signals are discussed in Appendix [A"l Two tempo- 
rally well separated incoming X-ray pulses interact with 
the molecule and excite two core-electrons at m and n 
(m ^ n) to the valence orbitals. The kjjj process is 
depicted by double-sided Fynmann diagram in Fig. (TT]) 

The corresponding correlation function expression is 
given by 



Sin(U, t 3 , t 2 , h) = 6{U - h)d{t 3 - t 2 )6(t 2 - h) 

[ (B m (i4)B„(t 3 )4('2)4(*i)) + (B n (U)B m (t 3 )Bl(ta)B^(ti))+U^t3 



E 



(5) 



The core exciton operators B and B^ are defined as 



Br, 



E^ 



Hi = E 



(6) 



The time evolution is in the Heisenberg representation, 



B{t) = e zH,t Bc- lH,t 



(7) 



where H' = Hq + He- £4 <=> £3 inside the brackets in 
Eq. ([5]) represents two terms obtained by an interchange 
of £3 and £4 in the first correlation functions and the 
expectation value in Eq. ([5]) is with respect to the N- 
electron ground state of the full (valence+core) system. 
For clarity we assume that the x-ray wavelength is short 
compared to the separation of sites n and m so that the 
phase variation of the field across the molecule may be 
neglected. This condition may not hold for hard x-rays. 
Equation ([5|) then need to be corrected by adding phase 
factors ea;p(ik.(r„ — r m ). The more general expressions 
are given in the appendix of [24| . 

A formally straightforward way to compute the sig- 
nal S is by expanding it in many-body states and their 
eigenvalues [20|. However computing the many-body 
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FIG. 1: The nonlinear signal in the direction = ki + 

k2 — k3 is depicted using the double sided Feynman diagrams, 
describing the evolution of the density matrix. The time in- 
creases from bottom to top. |o) is the ground state of the 
molecular system. |6) and |c) are the singly excited states 
while |e) is the doubly excited many-body state. 



states can become very demanding while efficient elec- 
tronic structure codes exist for valence excitations, these 
are much less developed for core-excitations. Here we use 
an alternative approach and express Sju in terms of the 
one and two-particle Green's functions. An algorithm is 
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presented for their computation which avoids the explicit 
computation of the many-body states. 

We are interested in the dynamics of two valence elec- 
trons in presence of the Coulomb potential generated by 
core-holes. The core- holes are treated as fixed parame- 
ters. 

We next introduce the one- and two-electron Green's 
functions [27112^, 



-i(Ta(t)c}(t')}m 

(- i ) 2 (r Cl (t) C ,(t')4(ti)4(t'i)) r 



(8) 



where (•)„, ((-) m)n ) denote the trace over N-electron 
ground state in the Coulomb potential of one (two) 
hole(s) at m(m,n). T is the time ordering operator; 
when acting on a product of time dependent operators, 
it rearranges them in increasing order of time from right 
to the left. The superscripts (m) and (m,n) represent 
parametric dependence on a single core-hole located at 
m and two core-holes at m and n, respectively. G\j can 
be expressed as, 

G^V) = -ie(t-t'){a\ci(t)c](t')\a) m (9) 

+ i9(t' -t)(a\c](t')ci(t)\a) m 

= -i6(t - t'y^-t'HalaUit - t')c)\a) m 

+ i6(t' -t)e- iE «( i -*')(a|ct[/(t' -t)a\a) m 

where a) is the valence many-body state with energy 
E a . The time evolution of operators in @ is given by 
a parametrized Hamiltonian, U(t) = e - lff ( m )*. H(m) is 
obtained from U)-© by tracing over core degrees of free- 
dom after replacing W imnj = 6 mn W-™^ and V mnm > n ' = 
Smn'Snm'V^ 11 ^ where l/( m >") is a parameter and 
depends on m parameterically. Note that W" 1 '™) = 



for m = n. Since the core excitations are stationary, the 
core operators are assumed to be time independent. The 
parametrized Hamiltonian thus obtained is hermitian. In 
writing the second equality sign we have further ignored 
the dependence of E a on the core-hole at to. Thus all 
dependence to core-hole is through the time evolution 
inside the trace, as marked by a subscript to. 

Similarly by parametrizing the core-hole, we can ex- 
press two-particle (two-time) Green's function as 

G\™f(t,t,t\t')= (10) 
(-i) 2 9(t - t'^t-^ialciCjUit - t')clc\\a) m ,n 
+ (-i) 2 9(t' - ^e^'-^ialclcjUit - t')c iCj \a) m , n 

where m =/= n. 

In order to express the signal in terms of the product of 
one- and two-electrons Green's functions, we need to re- 
duce our many-body space with (N+l) valence electrons. 
This is done by introducing the following projection op- 
erator in the (N + l)-valence/l-core-hole space 



|a)(a|<4 



(11) 



where \a) is a many-body state of the full (valence+core) 
system. This projection operator spans a subspace of 
the many-body space of N + 1-valence electrons (i.e. CI 
singles. The full space includes all possible multiples of 
excitons). This projection allows us to obtain a closed 
expression for the signal in terms of the one- and two- 
electron Green's functions, which can then be computed 
using many-body techniques. 

We next insert the projection operators (fTTj) into the 
expectation values in Eq. ([5]). This is our key approxi- 
mation which allows us to express the signal in terms of 
the product of Green's functions. Using Eq. ^ in (JSJ) 
and inserting the projection operators we obtain 



J 



x e 



Sin(t4, t 3 , t 2 , ti) — ^ ^ fJ-im^7nk^jn^nl&(t4 ~ ts)6(t 3 — t 2 )9(t2 — t\) 

m^n ijkl 

[(a| Ci 4C/(t 4 - t 3 )T mCj ciU(t 3 - t 2 )c m c\v n U(t 2 - h)c n c\\a) 
+ (a\c jC lU{U - h)V nCl c] n U{h - t 2 )c m clv n U(t 2 - h)c n cj\a) +t 4 ^t 3 

I 



The exact expression for the signal may be obtained by 
removing projection operators inside the brackets in Eq. 

Using (fTTj) . each term on the r.h.s. of Eq. (fT2"|) fac- 
torizes into a product of three correlation functions. By 
parametrizing the core variables, each correlation func- 



(12) 



tion can be approximated as 

(a|cJ„CjJ7(*4 - h)c\cm\a) w (a\ciU(U - h)c]\a) m (13) 

where U(ti — £3) on the r.h.s depends on core variables 
parametrically. 

Making use of Eqs. ©, (HHD and |T3]), Eq. (HJ) be- 
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comes 



Sni(t 4 , h, t 2 M) = 6{U ~ h)6{h - t 2 )9{t 2 - h) 

Mim Mrofe M jn f^nl Gj" (^2 1 £ 1 ) 

x [ Gg° (t 4 , t 3 )<#$ (t a , f a ) + G$ (t 4 , f a )<?## (*3,* 2 ) 

+ (*4<»t3) ■ 

(14) 

This is a closed expression for the signal in terms of the 
one- and two- electron Green's functions, corresponding 
expressions for other techniques Si at k/ = — ki +k2 + k3 



and Si i at k// = ki — k 2 + k 3 are given in Appendix [A] 

Note that the Green's functions in Eq. (|14p depend 
only on the ddifferences of their time arguments [see Eq. 
©]. Denoting the time delay between x-ray pulses as 
Ti = ti + i — ti, where i — 1,2,3, and taking the Fourier 
transform with respect to delay times, the signal can be 
expressed in the frequency domain as 



Sih(lu 3 , uj 2 , loi) 



we get 



oo />oo 



oo J —oo 



drx \ dr 2 I dr 3 c 



Siii{ti+t 2 +t 3 ,ti+t 2 ,ti,0), (15) 



SlIl(u>3,CJ2,Ul) 



m^n ijkl i'j' 



dLo[duj 2 duj 3 



LU 3 - W, 



(wi 



ie)(u> 2 



ie) 



LU3 



■ LJi — u n — ie 



(16) 



r 



where e is an infinitesimal number and 

/oo 
-oo 

is the Fourier transform of the Green's function. 
III. COMPUTING THE GREEN'S FUNCTIONS 



A self-consistent solution of the exact equations is highly 
demanding numerically. The GW approximation is gen- 
erally used in order to simplify these equations by ig- 
noring the vertex corrections and keeping only the first 
term in Eq. (|B31|) , the self-consistent calculation of the 
Green's functions and the self-energy is greatly reduced. 
We then have A^™] (t, t', h) = <5 ifc <%<5(i - ti)S(t - t') and 
Eqs. (|B29)l and ([Mil reduce to 



The one-electron Green's function, which appears in 
Eq. (fT4"|) , can be obtained from the self-consistent solu- 
tion of Dyson equation, 

G[f(t,t') = G%{t,t') 

+ G Q ijl {t^ x )tfj^ u ^)Gf;Ht'^nm 

where G° is the reference Green's function corresponding 
to non-interacting electron system and 

where Vq^ is the Coulomb potential due to the pres- 
ence of a core-hole at to. £( m ) is the time-dependent 
self-energy which comes from the electron-electron inter- 
actions and satisfy closed exact equations (|B28[) - (|B31[) . 



4- (mo = ^.tfMflW) (20) 

Si$(t,t') = V lkl] 8{t-t') + W vlkr S^, v {t' u t') 

x G^(t,t' 1 )G^)(t' 1 ,t) (21) 

where Vy/U — V%jki — Vjiki ■ Thus for one-electron Green's 
function we only need to solve Eq. ([T5)) self-consistently 
together with ([20]) and ([21]) . This is a standard approxi- 
mation which has been used extensively in studying opti- 
cal properties of different kinds of systems ranging from 
single atoms to semiconductor clusters and periodic crys- 
tals. The FEFF code[29( computes one-particle Green's 
functions by solving the GW equations, Eqs. ([T5|) - ([2"Tj) . 

The two-electron Green's function, Eq. ([8]), satisfies 
the equation (see Appendix iBj) 



+ Gp n \t,t' 2 )Gft n \t' 3 ,t')E%$^ (22) 
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where the exact interaction kernel, S, is given by Eq. 
(IB36[) . Note that one-electron Green's functions enter- 
ing in |22|) is calculated in the presence of two core- 
holes at m and n. These satisfy equations similar to 
Eqs. ([15]). P0|) and (f2"Tj) which can be obtained by sim- 
ply replacing superscripts (m) or (n) by (m, n) and us- 
ing E( m '™) = V^ m) + v£ n) + £( m >"). Recently Feng et 
al [30l ] have used nonresonant X-ray Raman scattering 
to study exciton spectroscopy on h — BN near boron 
if-edge by computing the two-particle Green's function. 
Within the GW approximation, the interaction kernel for 
the two-electron Green's function, Eq. (|B36|) . reduces to 
a simpler form, 

(*, t', t u t\) = xSgf V, t)S(t hW x 

(23) 



one-particle Green's function which includes the effect of 
the change in screened potential which goes beyond the 
usual GW approximation where variation in screening 
potential is generally ignored. The present formulation 
can be generalized straightforwardly to accomodate in- 
elastic interactions on the X-ray nonlinear signals by in- 
cluding their effects through self-energy for one and two 
particle Green's functions. 
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This is our final result for the scattering matrix be- 
tween two valence electrons in presence of two core-holes. 
The first term represents the contribution from screened 
coulomb interaction and the other two terms are induced 
by the change in screening of the electron in valance 
and hole in the core regions. Equations (f2"2"| and (j2"3"|) 
together with Eqs. dHJ), and (J3TJ) (after making 

above changes in superscripts for two core-holes) con- 
stitute a closed set of equations which can be solved 
self-consistently to obtain the two particle Green's func- 



tion. Since both S, 



(m.n) 



and G 



(m,n) 



come from the self- 



ijkl ik 

consistent calculation of the Dyson equation in presence 
of two-core holes, can be readily computed from 



Eq. (|23p . The numerical effort involved in computing the 
two-electron scattering-matrix is comparable to comput- 
ing the one-particle Green's function. 



IV. CONCLUSION 

We have derived closed expressions for the nonlinear 
X-ray spectra in terms of the many-body Green's func- 
tions. We expressed the k/// = ki + k2 — k3 signal 
in terms of one- and two-particle Green's functions. In 
the deep core Hamiltonian formulation, the slow core- 
hole dynamics is ignored as compared to the fast time 
evolution of the valence electron system, and core-holes 
are simply treated as parameters. The key approxima- 
tion was the introduction of a projection operator inside 
the multipoint correlation function of exciton operator 
which allows to express the signal in terms of one and 
two-particle Green's functions in a simple way. In or- 
der to compute these Green's functions, we have gener- 
alized the Hcdin's equations [l^] and derived a modified 
Bethe-Salpeter equation for two-particle Green's function 
in presence of two-core holes. A simple expression for the 
scattering matrix for two electrons is derived in terms of 



APPENDIX A: kj AND k /7 SIGNALS 



The k] — ki + k 2 + k 3 and k// = ki — k 2 + k 3 pro- 
cesses represent the interaction with two X-ray pulses and 
involve four Liouville space pathways which can be ex- 
pressed in terms of double sided Feynman diagrams (3l| 
shown in Figs, and ©. The corresponding expres- 
sions in terms of the correlation functions of the exciton 
variables analogous to ([5]) are given by 
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FIG. 2: Double sided Feynman diagram for the nonlinear 
signal in the direction kj = — ki + k2 + k3. Time increases 
from bottom to top. |5) is the ground state of the molecular 
system. \b) and |c) are the singly excited states while |e) is 
the doubly excited many-body state. 
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Similarly for kjj we obtain, 

m^n a ijkl i' j' 



x 6{U - t 3 )9(t 3 - t 2 )6(t 2 - t x ) 



^ " J {fc;li'f C ;i 



X (G$ (t 3 , ix)^ (t 4 , fe) + (*3, h)0Z$ (*4, *3))](A4) 

where / c; ij are the expansion coefficients in terms of the 
configuration interaction singles Q. 



APPENDIX B: DERIVATION OF EQS. (JH]) - ([22]) 



FIG. 3: Double sided Feynman diagrams for nonlinear signal 
in the direction Kn = ki — k2 + k3. 



Si(h,t 2 ,t 3 ,u) = ^e{u-t 3 )e{t 3 -t 2 )6{t2-tx) 

in n 

[(B n (ti)Bt(t 2 )B m (t 4 )Bt,(t 3 )> 
+ {B n {tx)Bl{t 3 )B m {U)Bl{t2)) 
+ (B n ( tl )Bl(U)Bl(t 3 )Bl(t 2 )) 
+ (B^ti)iM^(*3)i4(*2)>](Al) 



and 

Sn(t\,t 2 ,t 3 , ti) 



0(U-t 3 )9(t 3 -t 2 )9(t 2 -h) 

[{B n {t 2 )B\{t 3 )B m {n)Bl{t{)) 
(B n (U)Bi(t 3 )B m (t 2 )Bl(ti)) 
(B n (t 2 )B m (U)Bl{t 3 )Bl{tx)) 
(B n (t2)B m (U)Bi(t 3 )Bl(tx))](A2) 



where t^i = 1,2,3,4 is the interaction time of the four 
X-ray pulses, and operator B m is defined in Eq. 
Repeating the steps that transformed Eq. ([5]) to (fT4)) . 
we can express these signals in terms of the one and two- 
particle Green's functions as 

Si(h,t 2 ,t 3 , u) = J2J2 E E 

f-^im H>mj l^kn l^nl 

mn a ijkl i'j' 



x 9(U - t 3 )d(t 3 - t 2 )9(t 2 - h) 



We use a Hamiltonian parametrized in terms of the 
core-holes obtained by tracing out the core-hole degrees 
of freedom. 



i ijkl 



( m U„._ (B1) 



The first two terms represent the kinetic energy and the 
electron-electron interaction, </>™ in the last term rep- 
resents the potential due to the presence of a core-hole 
located at mth core orbital. It is obtained by approxi- 
mating the last term in Eq. ^ with, 



4r 5 



■WimmjCiCjiCmcln). (B2) 



We shall treat the valence electrons in the (average) field 
of two core-holes. We are interested in the dynamics of 
one and two electrons created by excitation from the core 
orbitals. We thus need to calculate 



Xwit^'Wi) = (Tc i (t)cJ(t , )c fc (ti)cJ(t / i)). (B3) 

which for t,t\ > t'i,t' represents the dynamics of two 
electrons created at times t\ and t' and destroyed at times 
t and t'\. Xijkl can be connected to a four-point Green's 
function <7^ n) as 

ftjfci(M',M'i) =iXi j kl{t,t\t 1 t\)+iG ij (t,t')G M {t\M m ) 

In order to obtain a closed equation for and fijS** , 
we add the following source term to the Hamiltonian, Eq. 

(ED. 



h g = Y ^ + E ^* c i 



(B5) 



2-*4) 



x G$\t 2 ,t 1 )G^{U,t 3 ) + t 2 &t 3 }-G%r(U,t 1 ) 



(«)t/ 



Here rj% and rj* are Grassmann variables which satisfy, 



(G$ (t 3) t a )0^ ia) + (i 3! t 3 )0£ft (t 4) t 3 ))](A3) 



_d__d_ 

% Or]* 



d d 
d-q* % ' 



(B6) 
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We further define a generalized one-particle Green's func- 
tion 



GijM := -——(Tc t (t)c](t')) v 
1(77,77*) 3 



(B7) 



where (•},, represents the trace with respect to the density 
matrix corresponding to the full Hamiltonian including 
the Grassmann terms and 

7(77,77*) = (Te- i ^ dT ^ Vi(T)c ^ T)+v * (T)c ^ T)) ) (B8) 

The one-particle Green's function G^' n ' is obtained sim- 
ply by setting rj = if' = in Eq. (|B7[) . Switching over to 
the interaction picture with respect to the Hamiltonian 
Hg, Eq. (|B7|) can be expressed as, 



Gij (t,t') 



T Cl (t)cUt')I( VlV *)) (B9) 



1(77,77*) 

Taking the derivative with respect to 7/ and 77* , we obtain 



dr ) l(t' 1 )dr 1k (t l ) 



= Gijkl(t,t',t 1 ,t'i) (BIO) 



M=o 



where {77} = is the short-hand notation for 77 = r/* = 0. 

Thus our strategy is to first compute a closed equation 
for the Green's function Gij and then using Eq. (|B10jl . 
we obtain a closed equation for the four particle Green's 
function, Q ijk i. 

The equation of motion for Gij is derived by taking 
time derivative of Eq. (|B7|> . 

Ag 

'dt 



7/(77, 77*)-G y (M') := S(t~t')5 l3 



T(- Ci (i)) C }(i'))(Bll) 



The time evolution of Ci(t) is governed by the total 
Hamiltonian H + Hq through the Heisenberg equation 

ij^(t) = e lCl (t) + cf%; n c f (t) + 2 V *,(t)ci(t)cl(t) 



2Tii'(t)ci'(t)ci(t) - Vi'ifk'cl(t)cf(t)c k '(t) - r\* 



(B12) 



wh ere V ijkl = {Vijkl ~ V jik i)/2. Substituting Eq. (|B12|) 
in (IBlip . we obtain, 



1/(77, rf) 



+ ia SjH + i^; n {t) 



G j , j (t,t') = 6 ij 5(t-t')+ir ) *(c](t')) TI 

iV iHj , k ,(Tcl,(t)c f (t)c k ,(t)ct(t% - 2i V *,(t)(Tc i (t)cl(t)c](t')} ri + 2i m ,{t){Tc v [t)c i {t)c]{t')) v . (B13) 

I 



We next connect the nonlinear terms in Eq. (|B13j) to the 
one-particle Green's function. By taking first and second 
derivatives of Gij with respect to 77^, we get 



5G k>j {t,f) 

+ 

S 2 G klj (t,t') = 



-(Tc r (t)c k ,(t)c)(t')) v 



(Tc k ,{t)c}{t')) v ( Cr {t))v 
i(Tc\,{t)c r {t)c k ,{t)c\tf)) n 
iG k , 3 {t, t')G jH >(t,t + ) 

1 s Vj ,(t) {cAt))v 



.5G k ,j(t,t'] 



s v m {Cj,{t)) " v 

Using Eqs. (|B14|1 in (|B13| we finally get 

G° fe ,(Mi)G^(*i,*') = 8 tJ 6(t-t') 

+ S^(i,ti)G^(ti,0(B15) 

where 



(B14) 



S(t - ^B16) 



with 



(m,n) 



(0 



2777*, (t) f(4(% + 
2ir) V {t) ((ci'ifyr, 
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S 



Svi'(t) 



+ (t>%' n \t)-i% ijlk ,G jH ,{t,t + ), (B17) 

Note that for {77} = the first two terms in (|B17|1 vanish 
and 7;™™ is the sum of Hartree and Coulomb potential 
due to two core-holes. The function Ey in (|B15|) is 



£*(Mi)G fc ,,(ti,i') = -ivt(t){cUt')) 



+ Vi 



i'ij'k' 



5 2 G k , 3 (t,t>) 



SG klj (t,t') 



:4(*)>. 



- SG k > j (t,f) / ... 



(B18) 



When {77} = 0, all terms, except the second term which 
is given by Eq. (|B10[) , vanish and E reduces to the self- 
energy E due to electron-electron interaction. 

In order to get a closed set of equations for the self- 
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energy S, we make use of the identity, 

G^(t,t 1 )lG fj (t 1 ,t') = 6 ij 6(t-t'). 



(B19) 



Differentiating Eq. l|BT9|) once with respect to %'(t), we 
can write 

= -S^G fe ^(M'). (B20) 

Using (|B20|) in (|B18j) , and doing some algebra we obtain 

Sij(t,i) = Wik'i> 1 j> 1 (t'it,t)Gk'i 1 Ci 1 ji' 1 j' 1 (ti,t'i 1 t'i) 
- ir ? *(t)(4(t' 1 ))^G^ 1 (t' 1 ,tO + A j (t,t') 



where 



Wijfcj(fi, £2, £3) = V»/yj 



S 2 v k i(h 



(B22) 



^(*3)^'(*a) 
-1 

y ' v = (B23) 



Svm (t 3 ) 



which for {77} = reduce to the screened Coulomb po- 
(B21) tential and the vertex function [l2|, HH, respectively, and 

I 



Aij(t,t') — iVi>ij>k>Gk>j 1 (t,t'i) 
- Vi'ij'k'Gkn^h) 
+ Vi'ij'k' 



SG^At'^t') . SG^At'^t') 
5 2 Gr)(t u t>) 5v, iri (t' 1 )8v ll3l {t' 2 ) 



Sr,*(t) 5r, f (t) 



+ Vi' Vj 1 



(B24) 



When {i]} = 0, Aij vanishes and only the first term on tively and making use of the identity (IB19|) . we obtain a 
the rhs of Eq. (IMT]) survives. closed equation for Cijki and Wijki 

Using Eqs. fB15j and (JB17J) in (|B23| and (B22} respec- 



&jkl(tl,t2,t 3 ) = —6ikSljS(tl — t 3 )S(tx - tz) 



(B25) 



— — Gi'j'^t i,t sjCj'J'^klit 3,t 4,,t 3 )Gj' 2 jAt i,t 2) 



5G l>r {t' u t'2) 

W ijk i(ti,h,t 3 ) = Q i jki( t ut2,t 3 ) (B26) 

— iVi' 1 kj' 1 l'Wiji> 1 j' 1 (t 3, t2,t 3 )Gj' 1 i 1 (t, t )£i 1 k 1 i'ij' 1 (t l,t 2, t' 3 )Gk 1 i> 1 (t' 2 , ti) 



where 



Qijkl{tl,t2,t 3 ) — iVi^kj'^Vi'ij'jGj'^^iA'l) 



SG~X(t' iA'2) , SGk^it^h) 



Sr)Ut 2 ) 8 Vj ,(t 3 ) 



5vUh)Sr]f(t 3 ) 



£jw + c.c. 



(B27) 
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On taking the limit {rj} = 0, Eqs. (|B15|) . (|B21j) reduce 
to a closed set of equations for the Green's functions and 
the self-energy 



(£,£') = G^ m> (t,t>) 

+ G^\tM)^,(hM)G ( ^(t 2l t) (B28) 

WO = s£? y (*i*,t)Gfc ) (ti,fa)Aja,(fa ) *i,*') 

(B29) 



0(m), 



where 5* and A are the screened Coulomb interaction and 
the vertex function, respectively, obtained from the re- 
duction of Eqs. {MS]) and (fB26f 



— 7-T— —G V j, (£ l,t 3 )A., ■, M {t 3 ,t 4,t 3 JG -, ,,(£ 4,t 2 J 



+ 



3' if 2^ 



3 23 



(B30) 



Sffi^i.ta.ts) = -iW^U^.ta)^ (B31) 

I 



The set of Eqs. (|B28|) - (|B31|) is exact. where 
We next derive a closed equation for the two-particle 

Green's function in the presence of a core-hole created 

at core orbital m. This can be readily generalized to the 

case of two core-holes. 

Equation (|B15|) can be written as 

Gij(t,t') — G°- (£,£') + Gjjtftjt' \)Y,j'k'(t'-i,t 'i)Gk'j(t' ' 2 ,£')- 



(B32) 



Differentiating Eq. (|B32|) . first with respect to rjkita) and 
then with respect to 77* (£4), taking the limit {rj} — and 
making use of Eq. (|B10[) . we get 



[<W<S(ti -£'2)-Gg, m) (£i,£'i)s5.r fe ) ,(£ / i,£ / 2)] 
+ G°j, m ' ) (t\,t'i)G k Tj (£'2, *2)3j-^ ife (£'i, £'2, £3, U) 



5 / <5S y (£ 1 ,£ 2 ) <5G f fc'//(t'i,t / 2)' 



_<fy*(*3) V^p(£'i,£' 2 ) ^(£4) 

JSij(£i,£2) (m) , / ,/ v 



J {7,}=0 

(B34) 



where in going from second to the third line we used the 
fact that, (SGij /8r)k){rj}=o — 0. 



Substituting Eq. (jB34)l in fl Bp) a nd using Eq. (|B28|) 
(B33) to replace the first line in Eq. (|B33|) . we get 



4«(*1,*2,*3,*4) = J Gi7 ) (£ 3 ,£ 2 )4 m) (^4) 



+ t'l)^™ (£'2, i 2)S l ( ™, ) , t , ife , i (£'1, £' 2 , £'3, *'4)£;jfcJ *2, £3, *4) 



(B35) 



r 



Eq uations (|B28| - ((B3T|l together with Eq. ||B35|) and 
(IB36[) constitute the set of equations (fTH|) - (|22|) used in 

,,„, / , , x S'E^(t'i,t'o) , , the main text. 

^V lfe ^'i,£' 2 ,£V 4 ) = -g ' \ (B36) 

dG ;' 1 fc' 1 (.r 3,£ 4j 



where kernel 
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